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Abstract 



We investigate the associated production of neutrahnos e~^e~ — > X1X2 accom- 
panied by the neutrahno leptonic decay X2 ~^ Xi'^~^^~^ taking into account 
initial beam polarization and production-decay spin correlations in the min- 
imal supersymmetric standard model with general CP phases but without 
generational mixing in the slepton sector. The stringent constraints from the 
electron EDM on the CP phases are also included in the discussion. Initial 
beam polarizations lead to three CP-even distributions and one CP-odd dis- 
tribution, which can be studied independently of the details of the neutrahno 
decays. We find that the production cross section and the branching frac- 
tions of the leptonic neutrahno decays are very sensitive to the CP phases. 
In addition, the production-decay spin correlations lead to several CP-even 
observables such as lepton invariant mass distribution, and lepton angular dis- 
tribution, and one interesting T-odd (CP-odd) triple product of the initial 
electron momentum and two final lepton momenta, the size of which might be 
large enough to be measured at the high-luminosity future electron-positron 
collider or can play a complementary role in constraining the CP phases with 
the EDM constraints. 
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I. INTRODUCTION 



The minimal supersymmetric standard model (MSSM) [Q] is a well-defined quantum 
theory of which the Lagrangian form is completely known, including the general R-parity 
preserving, soft supersymmetry (SUSY) breaking terms. The full MSSM Largarangian has 
124 parameters - 79 real parameters and 45 CP-violating complex phases . The number 
of parameters in the MSSM is very large compared to 19 in the standard model (SM). 
Therefore, many studies on possible direct and indirect SUSY effects have been made by 
making several assumptions and investigating the variation of a few parameters. Recently, it 
has, however, been shown [Q that limits on sparticle masses and couplings are very sensitive 
to the assumptions and need to be re-evaluated without making any of the simplifying 
assumptions that have been standard up to now. 

Despite the large number of phases in the model as a whole, only two CP-odd rephase- 
invariant phases, stemming from the chargino and neutralino mass matrices, take part in the 
chargino and neutralino production processes 0. In light of this aspect, phenomenological 
analyses with the complex parameter set are not much more difficult than those with the 
real parameter set in the chargino or neutralino systems. Incidentally, the CP phases are 
constrained indirectly by the electron or neutron electric dipole moment (EDM) and may be 
small, but the indirect constraints on its actual size depends strongly on the assumptions 
taken in those analyses. As a matter of fact, many recent works have shown that the 
constraints could be evaded without suppressing the CP phases of the theory. One option ^ 
is to make the first two generations of scalar fermions rather heavy so that one-loop EDM 
constraints are automatically evaded. This case can be naturally explained by the so-called 
effective SUSY models where de-couplings of the first and second generation sfermions are 
invoked to solve the SUSY FCNC and CP problems without spoiling naturalness. Another 
possibility is to arrange for partial cancellations among various contributions to the electron 
and neutron EDM's. Following the suggestions that the phases do not have to be suppressed, 
many important works on the effects due to the CP phases have been already reported; 
the effects are very significant in extracting the parameters in the SUSY Lagrangian from 
experimental data estimating dark matter densities and scattering cross sections and 
Higgs boson mass limits |T0|, CP violation in the B and K systems [|ll[], and so on. 

If the scale of the SUSY breaking is around 1 TeV as preferred by fine-tuning arguments 
in the Higgs sector, many sparticles are expected to be produced at future colliders such 
as the Fermilab Tevatron upgrade, the CERN Large Hadron Collider (LHC), or future 
e~^e~ colliders proposed by DESY, KEK and SLAC. Among the sparticles, non-colored 
supersymmetric particles such as neutralinos, charginos and sleptons are relatively light 
in most superymmetry theories. With R-parity invariance, charginos and neutralinos, the 
mixtures of the gauginos and higgsinos, are produced pairwise in e~^e~ collisions, either in 
diagonal or in mixed pairs. At LEP2 |T^, and potentially even in the first phase of e'^e~ 
linear colliders (see e.g. Ref. [0)) the chargino xf and the neutralinos Xi 2 be, for some 
time, the only chargino and neutralino states that can be studied experimentally in detail. 
Furthermore, as they are expected to be lighter than the gluino and in most scenarios lighter 
than the squarks and sleptons, those lighter chargino and neutralino states could be first 
observed in future experiments at e~^e~ colliders. On the other hand, the heavier chargino 
and neutralino states may require the second-phase e~^e~ linear colliders with a cm. energy 
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of about 1.5 TeV. 

In light of the previous generic arguments and aspects concerning the sparticle spec- 
trum, one of the most promising SUSY processes for investigating a wide region of the 
SUSY parameter space is the associated production |ll^ , p!5[| of neutrahnos in e^e~ colhsions: 
e^e^ XiX2- Although in general chargino production |T6| , p!7| is favored by larger cross 



sections, sizable cross sections for the neutralino process can be expected in certain regions 
of the parameter space. Moreover, it might be possible to discover SUSY by neutralino 
production if charginos are not accessible. 

If neutralinos as new particles are discovered, its clear identification can be enhanced 
by utilizing beam polarization and the complete investigation of the neutralino decays ||18|| . 
Polarized electron beams have been shown to play a critical role in disentangling SUSY pa- 
rameters in chargino-, neutralino- and sfermion-pair productions in e^e~ collisions. How- 
ever, those works have mainly considered longitudinal electron polarization. Recently, an 
intensive study to obtain high positron polarization has been made. In this light, we study 
the case where the polarization of both electron and positron beams can be freely manip- 
ulated, and we investigate if the highly polarized electron/positron beams can provide a 
powerful diagnostic tool for determining SUSY parameters in the associated neutralino-pair 
production. Angular distributions and angular correlations of the decay products as well as 
neutralino decay widths and branching ratios can give valuable additional information on 
their composition from gaugino and higgsino components. Certainly, one can infer the spin 
of the new particles from decay angular distributions with complete spin correlations of the 
decaying particle. Moreover, the identification of neutralinos can be very much solidified by 
ascertaining the Majorana character of the neutralinos [|T5|. This has been demonstrated 
to be possible by means of the energy distributions of the decay leptons if the neutralinos 
are produced in collisions of polarized e+e~ beams. The angular distributions of the decay 
products might, however, offer the possibility to prove the Majorana property although po- 
larized beams are not available. Furthermore, the angular distributions of the final leptons 



are suitable observables for studying CP violation in the MSSM [p!9[] . 

After the mixed pair x\ and X2 is produced, the second lightest neutralino X2 decays 
to the LSP and two fermions. Since leptons among fermions are most cleanly identified at 
high-performance detectors, one of the most promising modes for the associated production 
and sequential decays of neutralinos will be 

e+e" x\xl 

So, in the present note we give a comprehensive analysis through initial beam polarizations 
and spin correlations between production and decay to investigate the effects of the CP 
phases and the other SUSY parameters in the associated production and decay of the neu- 
tralinos in the MSSM with general CP phases but without generational mixing. In doing the 
analysis, it will be meaningful to include the stringent constraints from the electron EDM 
measurements on the CP phases Instead of performing full scans over the phases and real 
SUSY parameters, we take two typical scenarios to suppress the electron EDM constraints 
or to allow a large space of the CP phases while satisfying the electron EDM constraints. 
The choice of two scenarios will be made after a global study of the dependence of the elec- 
tron EDM on the relevant parameters in the MSSM. In each scenario, the neutralino mass 
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spectrum, three CP-even and one CP-odd obscrvablcs using initial beam polarizations, the 
ncutralino polarization vector and several distributions observable in the laboratory frame 
are presented and discussed. 

The organization of the paper is as follows. In Section II we describe the supersym- 
metric flavor-preserving mixing phenomena in a general parameterization scheme without 
generation mixing; selectron left-right mixing, chargino mixing and neutralino mixing, and 
identify the relevant CP phases. Section III is devoted to the discussion of the constraints 
by the electron EDM on the CP phases and Section IV to the associated production of 
neutrahnos with polarized beams and the introduction of useful CP-even and CP-odd ob- 
servables, which are extracted by controlling the initial beam polarization. In Section V 
we describe in detail the possible decay modes and branching ratios of the second-lightest 
neutralino. In Section VI, we explain how to obtain the fully spin-correlated distributions 
of the associated production and decays of the neutrahnos and study the impact of the CP 
phases on various angular correlations. Conclusions are given in Section VII. 



II. SUPERSYMMETRIC FLAVOR CONSERVING MIXINGS 

If flavor mixing among sleptons is neglected, the mass matrix of selectrons is given by 



z-2 



6 I /y-Ki f A _l_ I 1 * + -n ^ /KVi 2 _|_ 2 2 ^> r, O /-? ci2 I V / 



+ ml + ml cos 2/3 (s^ - 1/2) -nie (A* + n tan /3) 



* — me (^e + A**tan/3) + mg — m| cos2/3 



The first term of the diagonal elements is the soft scalar mass term evaluated at the weak 
scale, and the second is the mass squared of the the corresponding electron (dictated by 
SUSY), and the last comes from the T> term. The trihnear term Ag causing left-right mixing 
is due to the soft-breaking Yukawa-type interaction, and fi is the supersymmetric higgsino 
mass parameter describing the mixing of two Higgs doublets. And tan /? is the ratio V2/V1 of 
the vacuum expectation values of the two neutral Higgs fields which break the electroweak 
gauge symmetry. The selectron mass eigenstates can be obtained by diagonalizing the above 
mass matrix with a unitary matrix Ue such that UeAilUl — diag(m|^ , mf^ ) . We parameterize 
Uf. so that 

TT _ ( cos 6*6 — sin e"*''*''' \ /„s 
^~ l^sin^ee^-^^ cos^e y' ^ ' 

where —me{Ae + /i*tan/3) = \me{Ae + /x*tan/5)|e~*'^'= and we choose the range of 9e and (pe 
so that < < TT and -ti/2 < (/>e < ti/2. 

In supersymmetric theories, the spin-1/2 partners of the W bosons and the charged 
Higgs bosons, W"^ and H"^, mix to form chargino mass eigenstates Xi,2- The chargino mass 
matrix is given in the {W~ , H~) basis by 

A^<, = ( v/^"-"^"), (3) 

which is built up by the fundamental SUSY parameters; the SU(2) gaugino mass M2 as well 
as the higgsino mass parameter /i and the ratio tan (3. Since the chargino mass matrix A4c 
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is not symmetric, two different unitary matrices acting on tfie left- and right- chiral {W, H) 
states are needed to diagonalize tlie matrix: 



so that UrM-cUI = dia£r(m-±, m-±) with the ordering m-± < m-±. 

On the other hand, the neutral super symmetric fermionic partners of the B and 
gauge bosons, B and W^, can mix with the neutral supersymmetric fermionic partners of 
the Higgs bosons, and to form the mass eigenstates. Hence the physical states, x°, 
called neutralinos, are found by diagonalizing the 4x4 mass matrix 



M 



N 



f Ml -mzcpsw rrizspsw \ 

M2 mzC/sCw -mzspcw 

—mzcpsw rnzcpsw —/j, 

\ mzspsw -mzspcw -/J' / 



(5) 



where sp — sin/3, cp — cos (3, and sw,cw are the sine and cosine of the electroweak 
mixing angle 6w, respectively. The neutralino mass matrix A4n is a complex, sym- 
metric matrix so that it can be diagonalized by just one unitary matrix N such that 
N*AiNN^ = diagfm^o. m^o. m^o. m,^()) with the ordering m^o < m^o < m^o < m^o. 

Al' A2 ■ A3' X4 ^ Al X2 X3 A4 

In CP-noninvariant theories, the gaugino mass M2, Mi and the higgsino mass parameter 
II as well as the trilinear parameter can be complex. However, by reparametrization 
of the fields, M2 can be assumed real and positive without loss of generality since all other 
parameter choices are related to our choice by an appropriate R transformation. Taking into 
account our parameterization choice, the final set of phases considered in the discussion of the 
electron EDM and the neutralino production and decays includes two phases appearing in 
the chargino-neutralino sector $1, and one phase corresponding to the trilinear soft 
breaking parameter relevant in the electric dipole moment calculation as will be discussed 
in the following section: 

H = \ii\ e^*^ Ml = \Mi\ e'^\ A, = \A,\ e**^«. (6) 

We note that even though the off-diagonal elements of the selectron mass matrix are pro- 
portional to the small electron Yukawa coupling, the CP phases, in particular, play a 
crucial role in determining the size of the electron EDM because every SUSY contribution to 
the EDM requires a chirahty flip leading to dipole moments' proportionahty to the electron 
mass me- 



III. ELECTRIC DIPOLE MOMENT OF THE ELECTRON 

The electric dipole interaction of a spin-1/2 electron e with an electromagnetic field is 
described by an effective Lagrangian 

JO-edm = - 2^ ea'^'-f^e F^^. (7) 
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In theories with CP-violating interactions, the electric dipole moment de receives contribu- 
tions from one loop diagrams. In the MSSM, two diagrams contribute to the electron EDM 
in the mass eigenstate basis of all particles. They are shown in Fig. 1 (summation over all 
charginos and neutralinos in the loops is understood). 







Figure 1: Diagrams contributing to the electron EDM; (a) chargino- exchange contributions 
and (b) neutralino- exchange contributions. 



The Lagrangian describing the x^-e-z/g interactions without flavor mixing is 



(8) 



where Pl,r = (1 T75)/2, and the interaction Lagrangian describing the most general x°-e-e 
interactions are given in terms of mass eigenstates by 



V2s 



w 



5fPi + 5fPJX-ga + h.C 



,~,0 ~ 



where a = 1,2 and the couplings Bj and are given by 

Bj^ = V2Y^N*^ cos Oe + 2N*^ tan 9we''^' sin 9^ , 
B]^ = -{Nj2 + Nji tan 9w) cos 9^ + V2Y^Nj3e''^' sin 9^ 
Bf = -V2Y^N*^e-'^^ sin 9^ + 2N*^ tan 9w cos 9^ , 
= {Nj2 + Nji tan 9w)e-''^' + V2Y^Nj^ cos 9e , 



(9) 



(10) 



with the electron Yukawa coupling Y^ = rrie/ {\/2mw C/3) ~ 6.4 x 10~^/c/3. 

It is clear that the matrix elements of two unitary matrices Ul and Ur diagonalizing the 
chargino mass matrix Aic are functions of the phase but not of the phase $1. Using the 
chargino-electron-sneutrino interaction, we find that the chargino contribution to the EDM 
for the electron through the diagram shown in Fig. 1(a) is 



a 



2 I m-± 



m 



V 









Ulj2URjl 









(11) 
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where A{r) = 2(1 — r)~^[3 — r + 2 lnr(l — r)~^]. On the other hand, the neutrahno diagonal- 
ization matrix is a function of both $i and "l*^. Using the neutrahno-electron-selectron 
interaction, we find that the neutrahno contribution to the EDM of the electron through 
the diagram shown in Fig. 2(b) is 

4 ( m.-a ( \ 1 

(12) 

where Bir) = 2(1 — r)~^[l + r + 21nr(l — t)~^]. Our anlaytical expression for the electron 
EDM is consistent with that of Pokorski, Rosiek and Savoy of Ref. although there is a 
small difference in the neutrahno contribution between our result and that of Brhlik, Good 
and Kane P], of which the expression (AlO) must have a negative sign in the last term 
instead of a positive sign. These results are completely general except for flavor mixing and 
lead to the MSSM contribution to the electron EDM as the sum of two contributions: 

4 = df + df . (13) 

Of course, the Kobayashi-Maskawa CP phase in the SM can in principle contribute to the 
electron EDM, but it turns out to be effective only at three-loop level so that the contribution 
is too small to be measured. 

One of the important features of the SUSY contributions to the electron EDM is the 
fact that the EDM requires different chirality of the initial and final electrons. In the 
super symmetric diagrams this chirality flip can happen in two ways - either the exchanged 
selectrons change chirality via L-R mixing terms in the selectron mass matrix and couple to 
the gaugino component of the intermediate spin-1/2 particle, or the left- and right-handed 
selectrons/sneutrinos preserve their chirality and couple to the higgsino components of the 
charginos or neutralinos, respectively. As a result, all contributions are directly proportional 
to the mass of the external electron since both the L-R mixing selectron mass term and the 
Higgsino-electron-selectron (or sneutrino) coupling are proportional to the relevant electron 
Yukawa coupling Yg. Another consequence of the chirality flip is the explicit proportionality 
of the contributions to the mass of the intermediate spin-1/2 particle. 

Generally, the SUSY contributions to the electron EDM are determined by 7 real param- 
eters {tan/3, |Mi|,M2, mg^, mg^, \Ae\} and 3 CP phases {$i,$^,$Ae}- So, in order to 
understand the general features of the SUSY contribution effectively, it will be necessary to 
make some appropriate specifications without spoiling their qualitative aspects. We take a 
universal soft-breaking selectron mass mg for the left- and right-handed selectrons, because 
at any rate two states are not degenerate due to extra contributions to their masses. For six 
real SUSY parameters, we consider two typical scenarios where three CP phases are left as 
free parameters. In both scenarios, the gaugino mass unification condition is assumed only 
for the modulus of the gaugino mass parameters; \Mi\ = | tan^ 9w M2 ~ 0.5 M2 and the size 
of the trilinear parameter \Ae\ is set to 1 TeV through the paper. It is necessary to be careful 
in choosing the value of tan f3, according to which various physical quantities will be very 
different. The recent calculation by Chang, Keung and Pilaftsis pO| for the Barr-Zee-type 



two-loop contributions to the EDM's, the sbottom contributions are very much enhanced 
for a large tan f3 [^] so that the contributions cannot be simply neglected. Therefore, for 
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a large value of tan (3 we are forced to introduce more CP phases related with sparticles of 
the third generation in our analysis. Furthermore, as will be discussed in more detail in the 
section for the neutralino decays, for a large tan (3 we need to include stau left-right mixing 
as well as Higgs-exchange diagrams in evaluating different branching fractions of neutralino 
decays. On the contrary, the value of tan/9 < 2.5 has been already ruled out by null results 



in the Higgs search experiments at LEP II [^. Postponing the detailed analyses related 



with the tan j3 dependence of the EDM's, the associated production of neutralinos and the 
branching ratios of neutralino decays to our next work, we simply take tan/5 = 3 in the 
present analysis and treat the stau contributions on the same footing as the other slepton 
contributions. 

With these several specifications on the SUSY parameters, the electron EDM is deter- 
mined by three real parameters {M2,|/x|,me} and three remaining phases 
The first scenario 51, which is based on the so-called effective SUSY model, decouples se- 
lectrons by rendering them extremely heavy without violating the naturalness arguments, 
but taking M2 and relatively small: 

SI: M2 = 100 GeV, = 200 GeV, mg = lOTeV, (14) 

where 10 TeV for mg is taken because it is large enough to suppress the selectron contri- 
butions (almost) completely. As a result, the present electron EDM measurements pH] of 



\de\ < 4.3 X 10 ^^e ■ cm do not put any constraints on the CP phases. The second scenario 
52 takes a small universal soft-breaking selectron mass, but a large value of 

S2: M2 = 100 GeV, = 700 GeV, mg = 200 GeV. (15) 

In this scenario, we expect that some cancellations among the CP phases are needed to 
suppress the electron EDM. Of course, the degree of the cancellations depends on the values 
of the real SUSY parameters, especially the higgsino mass parameter The reason why 
we take a large |yu| of 700 GeV is to allow a relatively large region for the CP phases 
and $1 while scanning the phase from to 27r. The fact that the allowed region of two 
phases increases with has been pointed out in the work by Brhlik, Good and Kane 0. 

We display in Figure 2 (a) the allowed range of versus at 95% confidence level in 
the scenario S2 for other phases sampled randomly within their allowed ranges. The overall 
trend clearly shows that for larger values of it is much easier to satisfy the electron EDM 
limits and any value of larger than 650 GeV allows the full range of Figure 2(b) 
shows the allowed region at 95% confidence level for the phases $^ and $1 in the scenario S2 
with = 700 GeV. Note that near the region for $1 = tt the phase can take any value. 
In our next analysis on the correlated associated production and decay of neutralinos, which 
are mainly related with the CP phases $^ and $1 but not with the phase ^a^i Figure 2(b) 
will serve as the basic platform for all the contour plots for production cross sections, total 
cross sections of the correlated process, the branching ratios, several CP-even and CP-odd 
observables and an interesting CP-odd (T-odd) triple momentum product of the initial 
electron momentum and two lepton momenta from the leptonic decays of the neutralinos. 



IV. ASSOCIATED PRODUCTION OF NEUTRALINOS 
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A. Production helicity amplitudes 



Although we are mainly interested in one production process e+e" — > X2X1) '^^ discuss in 
this section the associated production of every combination of neutralino-pair e'^e" — > x°x° 
— 1-4] on a general footing. Note that the chirality mixing of scalar electrons are 
determined by the very small electron Yukawa coupling proportional to the electron mass 
much smaller than the collider cm. energy (500 GeV) under consideration by a factor of 
about 10®. Therefore, the selectron left-right chirality mixing can be safely neglected in the 
associated production of neutralinos so that the trilinear term does not play any role in 
the high energy process unlike the case for the electron EDM. In this approximation, the 
production process e^e~ —>■ Xi Xj is generated by the five mechanisms shown in Fig. 3: s- 
channel Z exchange, t-channel cl^/j exchanges, and -u-channel cl^^i exchanges. The transition 
matrix element, after an appropriate Fierz transformation of the eL,R exchange amplitudes 



(16) 



can be expressed in terms of four generalized bilinear charges, classified according to the 
chiralities a, P = L, R of the associated electron and neutralino currents 



LL 



LR 



RR 



+ 



Qrl — 



Dz , 
2 



2 ^2 y Vv 2 



— DuLdLij-i 
1 



H — Y'-^ij + P^tRgRiji 
D7 

--3-^13 ~ ^uRQRij, 



(17) 



with t-, and u-channel propagators and the couplings Zjj, guj and gRif 



m 



I + irrizTz ' 



D 



tL,R 



DuL,R — 



(18) 



with s — (pe +Pe)^) t — iPe ~ Px°)^ ^ — (pe — Pyp)^. And, the combinations Z^j, g^ij 
and guij of the neutralino diagonalization matrix elements Nij 



2ij - 2 



1 



9Rij 



4s^ 



(19) 
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satisfy the hermiticity relations reflecting the CP relations 

^ij ~ ^jii 9Lij = dLjv 9Rij = 9Rji-, (20) 

SO that, if the Z-boson width T z is neglected in the Z-boson propagator Dz-, the bilinear 
charges also satisfy the same relations = Q^" with t and u interchanged in the 
propagators. The relation is very useful in classifying CP-even and CP-odd observables in 
the following. 




Figure 3: Five mechanisms contributing to the production of neutralino pairs in e^e anni- 
hilation, e'^e~ XiX'j- 

All physical observables (which can be constructed and measured through the production 



process) are expressed in a simple form by 16 so-called quartic charges which contain 
important dynamical properties of the process and are expressed in terms of the bilinear 
charges Qa^- These quartic charges are classified according to their transformation proper- 
ties under parity as follows: 



(a) Eight P-even terms: 



Q 



Q 



4 

2 
1 



mni' + \q1l\' + mL\' + \q1r\' 



4 

il 

2 

il 

2 

n 

7^ 



\Qll\'^ + IQiflP \Qrl\'^ \Qlr\'^ 



^RR^RL "T ^LL^LR 
^RR^LR "T ^LL^RL 
^RR^LR "T ^LL^RL 



r)ij r)ij* 

^RR^LL 

r)ij Qij* 
^RL^LR 



(21) 
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(b) Eight P-odd terms: 



4 

In 

2 
1 



\Qrr\'^ + IQ^lP 



QrrQi 



RL 

IQ^RrI'^ + IQl/jT ~ IQrlI"^ 

^RR^RL h!LL^L 
^RR^LR ^LL^RL 



^LL^LR 



4 

il 

2 

In 

2 

ii 

2 
J 

I 



^RR^LR ^LL^RL 



Q%rQ''* 



RR^LL 

ij ^ij* 
RL^LR 



(22) 



We note that these 16 quartic charges comprise the most complete set for any fermion-pair 
production process in e^e~ colhsions when the electron mass is neglected. On the other 
hand, the quartic charges defined by an imaginary part of the bilinear-charge correlations 
might be nonvanishing only when there are complex CP-violating couplings or /and CP- 
preserving phases like recattering phases or finite widths of the intermediate particles. So, 
if there are no CP-preserving phases, non-vanishing values of these quartic charges signal 
CP violation in the given process. 

Defining the production angle with respect to the electron flight direction by 6, the 
helicity amplitudes can be determined from Eq. (0). Electron and positron helicities are 
opposite to each other in all exchange amplitudes, but the and helicities are less 
correlated due to the non-zero masses of the particles; amplitudes with equal neutralino 
helicities must vanish only oc m^o I \fs for asymptotic energies. Denoting the electron 

helicity by the first index, the ^ and x° helicities by the remaining two indices, the helicity 
amplitudes T(cr; Aj, \j) = 2na{(j; XiXj) are given by 



(+;++) = - 


Q%R^/i-vl + Q%LV 




sin 9 , 






(+;+-) = - 




-) + Q'U {I - ri. 




(1 + cos 9) , 


(+;-+) = + 


Q%R^ii-v+)ii-v- 


-) + Q^iL^/(l+^^ 




(1 - cos 9), 


(+;—) = + 


Q%RA/l'r^'- + Q%LV 




sin 9 , 






for the right-handed electron beam, and 










(-;++) = - 






sin 9 , 







(23) 
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(- 
(- 
(- 



+-) = + 
-+) = - 
— ) = + 



Q'1r^{^ - - V-) + Q'U{1 + ^+)(1 + ri-) 



[1 — cos 9) , 
^1 + cos 6) , 



LR 



sin 



(24) 



for the left-handed electron beam, where ri± = A^/^(l, /i^, ± {fi^ — yu|) with /i^ = m|o/s 

and A(x, y, z) = + y'^ + — 2xy — 2yz — 2zx. The explicit form of the production helicity 
amplitudes have been obtained by the so-called 2-component spinor technique of Ref. pSj . 



If the arguments are not specified, the notation A stands for A(l,/if,/i^) in the following. 



B. Neutralino mass spectrum and production cross section 

1. Neutralino masses 

Before investigating various dynamical distributions in the neutralino processes, it will be 
worthwhile to see the dependence of the neutralino masses and of the gaugino composition 
of the two light neutralino states on the CP phases in two scenarios SI and S2. Figure 4 
shows the mass spectrum of the neutralinos Xi,2 the plane of the CP phases $^ and $i in 
the two scenarios; (a) 51 (upper figures) and (b) 52 (lower figures). Except for the region 
around $^ = 0, 27r in 51, the second-lightest neutralino mass m^o is (almost) independent 
of $1 in both 51 and 52, while the lightest neutralino mass m^o exhibits a very strongly 
correlated dependence on the CP phases. Note that fn^o becomes maximal at non-trivial 
values of $^ and $i in 51 with = 200 GeV. This feature leads immediately to the 
conclusion that Tn^o is strongly affected by a small value of while Tn^o is essentially 
determined by the SU(2) gaugino mass M2. Combined with the electron EDM constraints 
shown as the shadowed region in Fig. 1, the mass m-j^o becomes smaller as the CP phases 
$1 and $^ approach the off-diagonal line on the plane, which implies that the mass is a 
function of the sum $^ + $1 of two CP phase to a very good approximation. 

Since the phase $1 is related with the gaugino part while the phase <I>^ with the higgsino 
part, the size of their contributions will be strongly dependent on the size of the gaugino 
(or higgsino) compositions of the neutralino states. So, we present in Figure 5 the gaugino 
compositions Xi and X2 of the lightest and second-lightest neutralinos Xi X2 defined by 





liVii 




|iVl2| 


X2 = 


|iV21 




|A^22| 



with respect to the phase $^ while the phase $1 is scanned over its full allowed range in 
the scenarios (a) 51 and (b) 52. As expected, Xi has larger gaugino composition than X2- 
Certainly, the gaugino composition in the scenario 52 is almost 100% due to the large value 
of |yu| = 700 GeV compared to the gaugino masses |Mi| and M2. For $^ around vr, the 
gaugino compositions are almost insensitive to the phase $1. On the contrary, in the region 
of $^ = 0, 27r, the gaugino compositions are very sensitive to the phase $1. This feature is 
partially responsible for the fact that in the scenario 51, the neutralino masses are strongly 
dependent on the phase $1 around $^ = 0, 27r. 
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2. Production cross section 



One of the important distributions in the production process e+e~ — > XiXj is the differen- 
tial cross section averaged over the initial beam polarizations. This unpolarized differential 
production cross section is given by taking the average/sum over the initial/final helicities: 

j^(.-e-^,?,»)^|gAV^EK-;A.A.)r. (26) 

Carrying out the sum, one finds the following expression for the differential cross section in 
terms of the scattering angle 6 and the quartic charges: 

X°X;) = ^AV^|[4 - (r^^ - ri_f + (r^+ + r^.f cos' e]Q'^ , 

+^^{i-vl){i-v'-)Q'^ + 4(7?+ + V-) cos eg^^j . (27) 

Figure 6 shows the dependence of the production cross section on the scattering angle 
and on the CP phases {$^, $1} in the scenario (a) SI and (b) S2 for a given cm. energy 
of 500 GeV. Several interesting features are noted: 

• Two distributions are forward-backward symmetric, which is due to the Major ana 
property of the neutralinos. This symmetry property can be traced back to the fact 
that the quartic charge Q*/ is directly proportional to cos© and the quartic charges 
Qh2 ^i^c forward-backward symmetric due to the Majorana relation \QaL Ri^)\ — 
IQai?,L(^ — 0)1 where a = L,R stands for the electron chirality. 

• The cross sections in the scenario S2 with small selectron masses are much larger 
in size that those in the scenario SI with very large selectron masses. This reflects 
the fact that the t- and -u-channel selectron exchanges become dominant for small 
selectron masses so that the production cross sections are very much enhanced. One 
additional crucial reason for the enhancement is that two neutralino states are more 
gaugino-dominated in the scenario S2 than in the scenario SI. 

• The cross sections due to the selectron exchanges are smaller in the forward-backward 
directions in contradiction with a naive expectation of forward-backward peaking phe- 
nomena due to t- and/or w-cxchanges. The reason is that while the Q^^ contribution is 
suppressed in the scenario 51 it is comparable in size with the other contributions with 
opposite sign in the scenario 52. Therefore, in the forward and backward directions, 
there exist a large cancellation among separate contributions. 

• The production cross section is more sensitive to the CP phases in the scenario 52 than 
in the scenario 51; the absolute production rate is much larger and the change due to 
different phases is larger as well in the scenario 52. So, we can expect in the scenario 
52 that the cross section itself can allow for a good determination of the phases. 

On the other hand the contours of the total production cross sections for the associated 
production of neutralinos e'^e" — > X1X2 displayed in Figure 7 on the plane of the phases 



da 



dcos© 



(e+e- 
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{$^, $1} for a given energy of 500 GeV and for two SUSY parameter sets of the scenarios 
(a) SI and (b) S2. The total production cross section is of the order of 1 fb in the scenario 
51 while it is of the order of 100 fb in the scenario 52. In the scenario SI, the cross section 
is large when the phase $1 is around tt and the phase is around tt/2 and 37r/2. However, 
the cross section in the scenario S2 increases as the phases and $1 approach the central 
point {tt, tt} along the diagonal as well as off-diagonal lines. 



C. Initial beam polarizations 

1. Spin-spin correlations 

At future e'^e~ colliders, it is expected that highly longitudinally polarized electron 
and/or positron beams are available. On the other hand, it is uncertain if high transversely 
polarized beams can be easily obtained unlike conventional e~^e~ circular colliders. Neverthe- 
less, it is interesting to investigate the effects of the longitTidinal and transverse polarizations 
of the initial beams for the determinations of the fundamental SUSY parameters. So, in this 
section, we introduce a general formalism to describe the polarization effects of the initial 
beams for any production process. Here, the extremely small electron and positron masses 
(me = 5.1 X 10^'^ GeV ) compared to the e~^e~ collision energy of the order of 100 GeV al- 
lows us to have a very much simplified formalism. Neglecting the electron mass renders the 
positron helicity opposite to the electron helicity in any theory preserving electronic chiral- 
ity as shown before. Let us consider the neutralino-pair production process e'^e~ — > x'iXp 
which is nothing but the process under consideration. For the sake of convenience we in- 
troduce a bracket notation for all the helicity amplitudes M^^-x^x = S„^^^{a : XiXj) which 
is guaranteed by the electronic chirality invariance and which enables us to obtain a simple 
form of the polarization-weighted squared matrix element as 



^ij = ^(1 - PlPl) E [K+ ^ ^^^.•>l' + K- ■■ M 



XiXj 

+^sin(a + «) E ■■ AiA,)(- : A.A,)*] , (28) 

XiXj 

where Pl{Pl) and Pt{Pt) denote the degree of longitudinal and transverse polarization of 
the electron (positron), and a{a) the direction of each transverse polarization with respect 
to a given reference plane, for which the scattering plane is chosen in most cases. The 
pictorial description of the azimuthal angles is given in Figure 8. We emphasize that the 
formalism can be applied to any e+e" collision process preserving electronic chirality with 
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an apporiate choice of reference frame to define the azimuthal angle parameters a and a. 



X 



Pt 




z 



Figure 8: Configuration of transverse polarization vectors in the center- of -mass frame. The 
azimuthal angles a and a denote the orientation of the polarization vectors with respect to 
a reference plane, for which the scattering plane is chosen. 

The longitudinal and transverse polarizations and the azimuthal angles for the transverse 
polarizations are related under CP transformations as follows: 

Pl " -Pl , 
Pt cos a < — > Pt cos a , 

Pt sin a < — > Pt sin a, (29) 

while the production helicity amplitudes are related under CP transformations as: 

{a : XiXj)ij < > {a : -Xj, -Xi)ji , (30) 

where the subscript ij means for the production of a particle Xi ^^id an anti-particle ac- 
cording to our spinor conventions. Denoting Ejj as the CP-conjugate polarization-correlated 
distribution of E^-, one can construct a CP-even distribution |(Sy + Sjj) and a CP-odd 
distribution |(Sjj — Sjj). If no CP-preserving phases are involved or any of them are neg- 
ligible in a given process, one can find that the CP-odd distribution is proportional to the 
last term of the distribution in Eq. (^), while the CP-even distribution is composed of the 
other three terms. 



2. Two additional CP-even observables 

Neglecting the Z-boson width {Tz ~ 2.5 GeV), which is very small compared to the 
collision energies under consideration, the distribution (|2^) provides us with three CP-even 
observables; one of them is the unpolarized part which has been discussed before and the 
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other two terms can be extracted by taking an appropriate polarization correlation. Longi- 
tudinal polarization yield a differential left-right (LR) asymmetry Alr 

^^« = ^E [K+;AA,)r-|(-;A.A,)r] , (31) 

XiXj 

with the normalization corresponding to the unpolarized part 



4 

The LR asymmetry can readily be expressed in terms of the quartic charges, 
Alr = {[4 - (r/+ - f + (7?+ + r}_y cos' e]g'f 

W(l-^+)(l-^-)^'2' + + V-) cos egf , (33) 

with, correspondingly, the expression for the normalization 

= [4 - (?7+ - 77_)2 + {rj+ + r]_f cos' e]Q\' 

+i^{l-nl) {l-ri'_)Qi^ + 4(77+ + ,7_) cos QQ'J . (34) 

It will be straightforward to extract the integrated left-right asymmetry experimentally with 
the expectation that highly longitudinally polarized beams are available at future e'^e" linear 
colliders. Certainly, the extraction efficiency depends linearly on the degree of electron and 
positron polarization obtainable at the e~^e~ collisions. 

The other transverse-polarization dependent term can be separated by allowing trans- 
verse polarization and setting longitudinal polarization to zero. Note that the transverse 
term is dependent on the sum of two azimuthal angles a and a, which must be sorted out by 
using a weight function \/2 cos(q; -|- a). This projection requires that the scattering plane is 
experimentally determined event by event. First of all, the neutralino masses m^o and m^o 
are expected to be measured with good precision through identifying the minimal and maxi- 
mal values for the lepton invariant mass in the leptonic decay X2 Xi^^^~- The determined 
masses and the four-momentum of two final leptons enable us to determine only the polar 
angle between the X2 flight direction and the Xi flight direction in the laboratory frame. 
Certainly, if the cm. energy is so large that neutralino masses are negligible, the neutralino 
direction can be identified with the direction of the two-lepton momentum. However, for 
a moderate cm. energy, it is not possible to completely determine the X2 scattering angle. 
Nevertheless, this distribution will affect the final two-lepton distribution partially so that 
it is not useless to investigate the dependence of the observable on the SUSY parameters. 
The CP-even observable Vt obtained through the angular projection procedure is given by 

= E ^ [( + ; ^^^^r] = © ■ (35) 

V XiXj 

The upper figures in Figure 9 exhibit the LR asymmetries A^ji and the lower ones the CP- 
even observables Vt for the production of the associated pair X2 Xi ^ ^ function of 
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the scattering angle at a cm. energy of 500 GeV in the scenarios (a) 51 and (b) 52 
for five combinations of the CP phases {$^,$i}. We note that these two observables are 
also more sensitive to the CP phases in the scenario S2 than in the scenario SI. However, 
the sensitivity of the LR asymmetry A^^ to the CP phases is not strong, so that the 
asymmetries are not very useful in determining the phases. On the other hand, as discussed 
before, the CP-even observable, which is more sensitive to the phases, is not easy to extract 
experimentally. Therefore, we may conclude that these two observables are not so powerful 
in determining the phases in both scenarios, while satisfying the constraints from the electron 
EDM measurements. 



3. One CP~odd observable 

CP violation arises in the existence of nontrivial complex couplings in the Lagrangian. 
In the associated production of the neutralinos CP violation is reflected in the complex 
production amplitudes. First of all, we find that in every diagonal production of neutralinos 
the production amplitude is purely real with the Z-boson width neglected, and it leads to 
no CP-violation. 

Like the CP-even observable Vt-, the only T-odd term Vn-, which is CP-odd in the absence 
of any CP-even rescattering phases like the Z-boson width, can be separated by allowing 
transverse polarization and setting longitudinal polarization to zero with A/2sin(a -|- a) as 
a projection angular function. As a result, one can obtain a CP-odd observable Vn as 



V 



1 



N 



.2v^A^sin2e. 



(36) 



One can check with the definition of the quartic charge Q'g that if the Z-boson width is 
neglected, the T-odd observable Vn niay be non-zero only for i ^ j since for i = j all the 
bilinear charges are real. 

In CP-noninvariant theories the quartic charge Q'g^, which is non-vanishing for i 7^ j, 
can be expressed in terms of two Jarlskog-type CP-odd rephasing invariants [2^ of the 
diagonalization matrix N . In order to elaborate on this point further, we present the explicit 
form of the quartic charge; assuming a real Z-boson propagator the quartic charge Q'g is 
given by 



Q 



nj _ 

6 - r, 4 4 



Sw{DtL - DuL)T{Zijgl^., 



l/2){D,n-D^n)l{Z,^g*) 



(37) 



Two combinations of the couplings, 2{Zijg1^j) and 2{Zijg^^j), are rephase- invariant. Using 
the expressions for Zij, guj and gmj, we can rewrite the two combinations as 



A^ijglij) 



2r2 



i{N,,N;,N*,N,r) -i{n,,n;^n:,n,, 
iiN,sN;,N':,N',2) - i{n,,n*,n;^n' 



(38) 



17 



where the primed matrix elements N--^ and are related with the diagonalization matrix 
elements and through 

= cw Nil + sw Ni2, A^2 = -sw + cw Ni2 . (39) 
Prom the expressions (37) and (38), we can draw the following consequences: 

• Both I{Zi.jg^.-.j) and T{Z,ijg*j^-) require the existence of gaugino and higgsino compo- 
nents and a different magnitude of two higgsino components of the x° and states. 
This latter requirement means that tan /3 should be different from unity. 

• The distribution is forward-backward asymmetric, because the angular dependence is 
determined by the difference DtL,R — D^l.r- 

• Due to the large suppression in the t- and -u-channel selectron exchanges, the T-odd 
asymmetry is very small in the scenario SI. Moreover, the asymmetry Vn is very 
small in the scenario S2 as well. This suppression is because the observable requires a 
sizable mixing between gaugino and higgsino states, but the mixing is very small due 
to the large value of |//| compared to the gaugino masses M2 and \Mi\. 

As a whole, these features lead to the conclusion that the T-odd observable Vn is not useful 
in measuring the CP phases directly in both scenarios SI and S2 suggested by the analysis 
for the electron EDM constraints. 

4- Neutralino polarization vector 

Neutralinos are spin-1/2 particles and their polarization can be measured through their 
decays. Before we investigate the possible neutralino decays in detail, in this section we 
study chargino polarization directly in the production process e+e" X2X1 with unpolarized 
initial beams. The polarization vector V"^ = {V^f^ ,Vqi ,V]^) of the produced neutralino is 
defined in the rest frame in which the axis z\\L is in the flight direction of x^, x\\T rotated 
counter-clockwise in the production plane, and y — z x x\\N of the decaying neutralino x°. 
Accordingly, the component Vi denotes the component parallel to the x° flight direction 
in the cm. frame, the transverse component in the production plane, and V^j^ the 
component normal to the production plane. These three polarization components can be 
expressed by helicity amplitudes in the following way: 

I ++)!' + ! +-)!'-! -+) I' - 1 {'^■,—)\'}/^f, 

(T=zh 

^t' = j:[{a;++){a;-+r + {a;—){a;+-r]]/Af, 

n = \r{ E K'^;— - (<t;++)(<7;-+)*] ]/n. m 

The longitudinal, transverse and normal components of the Xi* polarization vector can be 
easily obtained from the production helicity amplitudes. Expressed in terms of the quartic 
charges, they read: 
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Vi = 4|2(1 - /x^ - cos QQ'i' + cos QQ'^' 

+ Ai/2[1 + cos^ e - (/i^ - sin^ e]Q3^^j/Ar, 

P^^- = _8| [(1 - ^2 ^ ^2)g.y _^ yl/2Qnj q]^. ^ ^ ^2 _ ^2)^^.g^.jj q/^ ^ 

V% = ^\^'^^iJsinQQ'i/M , (41) 

where the reduced masses = m'i±/s. The longitudinal and transverse components are 

P-odd and CP-even, and the normal component is P-even and CP-odd. 

The normal polarization component can only be generated by complex production am- 
plitudes. Non-zero phases are present in the fundamental SUSY parameters if CP is broken 
in the supersymmetric interaction. Also the non-zero width of the Z boson and loop cor- 
rections generate non-trivial phases; however, the width effect is negligible for high energies 
as mentioned before, and the effects due to radiative corrections are small as well. So, 
the normal component is effectively generated by the complex SUSY couplings. As the 
selectron-exchange contributions can be ignored in the scenario 51, the CP-odd quartic 
charge simplifies to 

= (42) 



Since the presently measured value of s^r is 0.2315 '^^ry close to 0.25, the quartic charge 
Q4 is extremely suppressed in the scenario SI. However, in the scenario S2, the quartic 
charge can be relatively large without such a big suppression as shown in Figure 10. In this 
case, the normal polarization is of the order of 10 %, which is really sizable for non-trivial 
CP phases and very sensitive to the CP phases {$^, $1}. So, it is expected to give stringent 
constraints on the phases These strong constraints will be explicitly demonstrated in 
the following. 



V. NEUTRALINO DECAYS 

A. Decay density matrix 

Assuming the lightest neutralino Xi to be the lightest supersymmetric particle (LSP), 
several mechanisms contribute to the leptonic decays of the neutralino Xi > 2): 

X°te)-X?(go)+r(g)+£+(g) 

In particular, the leptonic three body decay of the second lightest neutralino, X2 ~^ Xi^~^^~ ■> 
is known to be very important because the end point of the lepton invariant mass distribution 
gives us direct information on the mass difference between X2 and Xi? which provides us with 
a stringent constraints on MSSM parameters. 

Although we will take into account only electrons and muons for the final state leptons, 
let us make some comments on the other possible leptonic decay of the second lightest 
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neutralino, X2 ~^ Xi'^~^'^~ IMIl- Since r is the heaviest lepton with a much larger mass (1.777 
GeV) than the other leptons and it couples with Higgs bosons with the strength proportional 
to tan f3, the branching fraction of this leptonic decay mode can be very different depending 
on the value of tan (3 and the Higgs mass spectrum. Actually, the mode is known to be very 
much enhanced due to the Higgs exchanges at large tan j3 [^. Furthermore, the polarization 
of r can be observed through the decay distributions, which strongly depend on the parent 
r polarization 



Due to the missing tau neutrinos, one would not be able to measure the invariant mass of 
two tau leptons experimentally. Nevertheless the r polarization or the invariant mass of the 
two r jets might be seen in future collider experiments. We note that in r ^ p or ai decays, 
the final vector meson carries a substantial part of the parent r momentum, therefore the 
smearing of the distribution is less severe than for decays into ir"^, fi and e. Let us give 
several comments on the tau decay mode. First of all, for f, the effects of the Yukawa cou- 
plings and slepton left-right mixing could be very important for large tan/3. Their leading 
contribution flips the chirality of the r lepton p9|. Secondly, for the three body decays, 
studying the correlation of two tau decay distributions would reveal the helicity flipping and 
conserving contributions separately. Thirdly, staus could be lighter than the other sleptons 
for various reseaons. The running of stau soft SUSY breaking masses from the Planck scale 
and stau left-right mixing could enhance decays into x'^t'^t'. Experimental consequences 
of such scenarios have recently been widely discussed. Also models with lighter third gener- 
ation sparticles have been naturally constructed without causing the flavor changing neutral 
current problem. The three body decay branching ratio and the decay distribution might 
be different from those for leptons in the first two generations. Since the study of the decay 
X2 ~^ Xi'^~^'^~ ill addition to the other leptonic modes could be an important handle to 
identify such models, we plan to present a detailed investigation about all these interesting 
features in the near future. 




Figure 11: Neutralino decay mechanisms; the exchange of the neutral Higgs boson is neglected 
because of the tiny electron Yukawa coupling. 

The diagrams contributing to the process Xi ~^ with i = e, fi are shown in Fig. 11 

for the decay into lepton pairs. Here, the exchange of the neutral Higgs boson [replacing 
the Z boson] are neglected since the couplings to the light first and second generation SM 
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leptons are very small. In this case, all the components of the decay matrix elements are 
of the left /right current x current form which, after a simple Fierz transformation, may be 
written for lepton final states as 



(43) 



with the generalized bilinear charges for the decay amplitudes: 

D' , o 1. 



^2 „2 y'^W 



Dlr 
Drr 



D'z 



Zli + A'ii 9Rli ) 



n' 1 



D' 



(44) 



where the chiralities a/f3 stand for the Xi/^~ chiralities and s'—,t'—, and m'— channel propa- 
gators and the couplings Zij, guj, and qru are given in the section for the production helicity 
amphtudes. The Mandelstam variables s', t', u' are defined in terms of the 4-momenta of 
X?,^~ and respectively, as 



s ={q + qf , t' = {qo + qf , u = (go + qf ■ 
The decay distribution of a neutralino with polarization vector n'^ is 



(45) 



P|2(n) = -4{t' - mlo){t' - mlo){Di - D^) - A{u' - ml,){u' - m^o)(Z)i + D 



,2 



-Sm^om^os'Do 



(n ■ q) m^o(m|() — u'){D[ + D'^) + m^o(m|() — t') 



D' 



+8{n ■ q) [-m^o(m|o - t'){D[ - D'^) + m^o{m\o - u')D: 
+16171^0 (q^nqq) D 4 , 



(46) 



where is the Xi spin 4-vector and {qinqq2) = ^fivpaq'i'nyq'^q" ■ Here, the quartic charges 
{DitoD^} and {D[toD'^ for the neutralino decays are defined by 



D2 



\D 



RR\ 



+ \DllV+\DrlV+\D_ 



'lr\ 



n [DrrDI^ + DllD*^^] 



\D 



LL\ 



RR\ 



-\DRLf-\D 



D, = -I[DrrDIj, + DllD},l] , 



\D 
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D' 



\D 



RR\ 



+ 1^ 



LR\ 



\D 



RL\ 



\D 



LL\ 



(47) 



If needed, the polarization of the final neutralino can be incorporated in a straightforward 
manner although the decay distribution will be more complicated in its form. 

For the subsequent discussion of the angular correlations between two neutralinos, it is 
convenient to determine the decay spin density matrix p\\i ~ V\Dy. In general, the decay 
amplitude for a spin-1/2 particle and its complex conjugate can be expressed as 

I?(A) = r uiq, A), V*{\') = uiq, X') f , (48) 

with the general spinor structure F and F = 7°F^ Then we use the general formalism to 
calculate the decay density matrix involving a particle with four momentum q and mass m 
by introducing three space-like four vectors (a = 1,2,3) which together with q/m 
form an orthonormal set: 



n 



9 



ab 



9ab n^nl 



9^u ■ 



(49) 



where g^''' = diag(l, -1, -1, -1) and g""^ = diag(l, -1, -1, -1) with a, 6 = {0 - 4}. A 
convenient choice for the explicit form of n°' is in a coordinate system where the direction 
of the three-momentum of the particle is g = (sin^, 0, cos^) lying on the x-z plane: 



n 



(0, cos6', 0, — sin 9) , 



n 



(0,0,1,0) 



n 



-M^Eq) . 



m 



(50) 



Then in this reference frame, n^'^'^ describe transverse, normal and longitudinal polarization 
of the particle. 

With the four-dimensional basis of normal four-vectors n^, n^, n^}, we can derive 
the so-called Bouchiat-Michel formula |^| 



u{q, X)u{q, A') = - [6xy + 75 J^'r^^x] 



m] 



which can be used to compute the squared, normalized decay density matrix p\\ 

V{\)V*{\')_ _ 1 
" ~ 2 



Paa' 



5xx' + -Y^l>x 



(51) 



(52) 



Tx\D{W 

where t"" (a = 1,2,3) are the Pauli matrices and the four kinematic functions X and Y"" 
(a = 1,2, 3) 

X = -8{t' - mln){t' - m\){Di - D^) - %{u' - m\){u' - m\){Di + D^) 



= -16(n" ■ q) 



m^o(m-o 



u'){D[ + D'^) + m^o(m|o - t')D'^ 



(53) 



+16(n'' ■ q) [-m^o(m|o - t'){D[ - D'.^) + m^o(m|o - u')D'^ 
+32m^o{qin°'qq)D4 , 

with n"- (a = 1,2,3) three vectors forming the polarization basis for the decaying spin-1/2 
particle. 
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B. Branching ratios 



Since the reconstruction of the neutrahno-pair production depends on the efficient use of 
the neutrahno decay modes, it is necessary to estimate the branching fraction of each decay 
mode. Since we assume that the hghtest neutrahno is the LSP and we are interested only in 
the decay of the second-hghtest neutrahno we can classify the decay modes as follows: 

^ ir,r^i>*,qr - xTr,xlqq- (54) 

Besides, if the mass m-± is smaller than the neutrahno mass m^o, the lightest chargino xt 
can enter the neutrahno decay chain via x'^ ~^ xt"^^* i xf H^* ■ Concerning the neutrahno 
decays, there are several aspects worthwhile to be commented on: 

• For the first and second generation fermions, the Higgs-exchange diagrams are sup- 
pressed unless tan (3 is very large. 

• The experimental bounds on the Higgs particles are very stringent so that the two- 
body decays x'^ Hxi and X2 ~^ -f^^xf are expected to be not available or at least 
strongly suppressed. 

• The lightest chargino and the second-lightest neutralinos are almost degenerate in the 
gaugino-dominated parameter space so that the charged decays such as X2 ~^ xt^^'^ 
will be highly suppressed. 

Nevertheless, we calculate the Icptonic branching fractions B{x2 ~^ Xi^^^ ) fully incorpo- 
rating all the possible decay modes of the neutrahno X2 while neglecting the Higgs-exchange 
contributions for a small tan/3 = 3. In our numerical analysis, we assume 200 GeV for a 
common soft-breaking slepton mass and 500 GeV for a common soft-breaking squrk mass 
in the scenario S2 while we take 10 TeV for all the soft-breaking sfermion masses in the 
scenario SI. Figure 12 shows B{x2 x"^"*"^") fo^^ £ = e or /i in the scenarios (a) SI and (b) 
52. We find that the branching fractions are very sensitive to $1 only around $^ = 0, 27r in 
the scenario SI, while it depends very strongly on $1 and on (almost) the whole space 
of the phases in the scenario S2. Furthermore, the branching fraction B{x2 ~^ Xi^"*"^") is 
very much enhanced in the scenario S2 because the slepton-exchange contributions due to 
mainly the gaugino components of the neutralinos become dominant due to the small slep- 
ton masses while the large value of suppressing the higgsino components. Consequently, 
the branching ratio of the leptonic decay of the second lightest neutrahno X2 Xi^^^^ is 
very sensitive to the values of the underlying parameters, in particular, the CP phases. 

VI. SPIN AND ANGULAR CORRELATIONS 

A. Correlations between production and decay 

In this section, we provide a general formalism to describe the spin correlations between 
production and decay for the process e'^e" — > x°x° followed by the sequential leptonic decay 
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Xi ~^ Xi^~^~^- Fo^^ the sake of convenience, wc do not consider the initial beam polariza- 
tion, which can however be easily implemented. Formally, the spin-correlated distribution 
is obtained by taking the following sum over the helicity indices of the intermediate neu- 
tralino state x° by folding the decay density matrix and the production matrix formed with 
production helicity amplitudes: 



-•unp 



X X X . 



(55) 



where the functions of the scattering angle O are given in terms of the production helcity 
amplitudes by 



-•unp 

V 
V 
V 



7 E I ('^; ++)!' + 1 ('^; +-)!' + ! -+)!' + 1 — > 

\ E [l ++) r + 1 +-) r - 1 ('^; -+> ? - 1 — ) 

\ E 7^{(a; -+)(<7; ++)* + {a- —){a; +-)*}, 

I E ^{(^; -+)(^; ++)* + — )(^; +-)*}■ 



(56) 



Notice that the above combinations are directly related with the polarization vector of the 
neutralino X2 as follows 



V 



V 



AT 



V 



unp ^unp ^unp 



(57) 



Combining production and decay, we obtain the fully-correlated 6-fold differential cross 
section 



da 
d$ 



where 



Sunp B(x° ^ X?r £+) 1 + P,V2 + P.V't' + PyK 



p -Yi. p - Yi p -Yi 

" X' ' X' ^ X 



(58) 



(59) 



with the phase space volume element d$ = dcos0da;ida;2dcos^^id0id0i2. The angular 
variable 9i is the polar angle of the l~ in the x° rest frame with respect to the original flight 
direction in the laboratory frame, and 0i the corresponding azimuthal angle with respect to 
the production plane, and 0i2 is the relative azimuthal angle of along the direction 
with respect to the production plane. On the other hand, the opening angle 612 between 
the ^~ and is fixed once the lepton energies are known. The pictorial description of the 
kinematical variables for the decay X2 ~^ Xi^~^^~ is presented in Fig. 13. The dimensionless 
parameters Xi and X2 denote the lepton energy fractions 



A2 



Ei+ 



A.2 



X2: 



(60) 
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with respect to the neutrahno mass m^o divided by a factor of two. The kinematically- 
allowed range for the variables is determined by the kinematic conditions 

< e < vr, 

< ^1 < vr , < 01 < 27r , < 0i2 < 27r , 

< 2;i,2 < 1 - r2i , (1 - - 0:2) > r2i , Xi + 0:2 > 1 - r2i , (61) 

where r2i = m-o/m?o, and the masses of the final-state leptons are neglected. 




Figure 13: Configuration of momenta in the e~^e cm. frame for the associated production 
of neutralinos and in the rest frame of the decaying neutralino X2- 

B. Total cross section of the correlated process 

The total cross section for the correlated process e~^e~ X1X2 ~^ is given 

by integrating the differential cross section (|58|) over the full range of the 6-dimensional 
phase space. More simply, the cross section is given as the multiplication of the integrated 
production cross section a{e~^e~ X2X1) the branching fraction i3(x2 ^ Xi^~^^~) 
because the correlation effects are washed out after integration over the phase space volume. 

The total cross section is displayed in Figure 14 as the contour plots on the plane of 
two CP phases {$^,$1} in two scenarios (a) SI and (b) S2. First of all, we note that 
the total cross section in the scenario 51 is very small. By definition, the scenario 51 has 
extremely large selectron masses so that only the Z-exchange diagram involving only the 
higgsino composition of the neutralinos contributes. However, the higgsino composition of 
the neutralinos is at most 40% and in most cases 20% for while it is at most 20% and 
in most cases 10% for Xi- Therefore, it is naturally expected to have a strongly-suppressed 
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cross section in the scenario. On the other hand, the total cross section is very much 
enhanced in the scenario S2, mainly due to large gaugino composition and small selectron 
masses which contribute to the cross section through the t- and w-channel exchanges. So, 
we can conclude that it is possible to have a relatively large cross section if the higgsino 
composition of the neutralinos is large or the slepton masses are small [^]. 

Quantitatively the scenario SI will present a few events of the neutralino process for a 
integrated luminosity of the order of 100 fb~^ while the scenario S2 give a few thousand 
events for the same integrated luminosity. So, a good precision measurement of the relevant 
SUSY parameters might be performed in the scenario S2 at future high luminosity e^e~ 
collider experiments such as TESLA, while it might be difficult in the scenario SI. As can 
be read from the figures, the cross section increases as the phase $i approaches vr in both 
scenarios. However, the dependence of the cross section on the phase is very different in 
two scenarios. In the scenario S2, the cross section increases monotonically as <l>^ approaches 
TT, but it becomes maximal at a non-trivial value of between (27r) and vr in the scenario 
SI. 

Compared with the dependence of the neutralino masses on the CP phases displayed in 
Figure 4, the cross section can be larger for larger neutralino masses and vice versa for quite 
large region of the CP phases. This implies that within the range allowed in the scenarios 
the signals with larger masses can have more possibility of being detected while those with 
smaller masses may escape detection. In this sense, future high luminosity experiments can 
give constraints on the CP phases simply by putting the upper limits on the event rates. 



C. Dilepton invariant mass distributions 

The invariant mass of two final-state leptons m«, which is nothing but the square root 



of the Mandelstam variable, v s', 

mu = Vs' = m|o {xi + X2 - I + r2i) , (62) 

is a Lorentz-invariant kinematical variable so that it is easy to reconstruct by measuring 
the energies of two final-state leptons Furthermore, the distribution for the invariant 
mass mil is independent of the specific production process for the decaying neutralino. This 
factorization is due to the fact that the invariant mass does not involve any angular variables 
describing the decays so that the polarization of the decaying neutralino is not effective. 

Figure 15 shows the two-lepton invariant mass distribution in the scenarios (a) SI 
and (b) S2. This distribution must reflect the two-lepton invariant mass distribution of 
the neutralino decay X2 ~^ Xi^^^~ itself multiplied by the total production cross section 
o"(e"*"e~ — > X2Xi)- As shown before, the leptonic branching ratio is very sensitive to the 
values of the underlying SUSY parameters. We flnd that the distribution of the invariant 
mass mil of the flnal-state two leptons is sensitive to the the CP phases as well. First of 
all, the end point of the maximal invariant mass is strongly dependent on the CP phases 
so that after all the real parameters are determined, the measurement of the end point will 
provide us with a very good handle to determine the CP phases. As one can notice from 
Figure 15, the sensitivity is larger in the scenario SI with a small parameter, which is 
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more comparable to the value of M2 than in the scenario 52. It clearly implies that the 
effect of the CP phases is enhanced for comparable gaugino and higgsino parameters. 

In passing, we note that since it is independent of the production mechanism, the lepton 
invariant mass distribution can be an important tool for studying supersymmetric models 
even at hadron colliders because of its clean signature. This point has been in detail explored 
by Nojiri and Yamada |^ by investigating the parameter dependence of the distribution of 
the three body decay xl ^ X?^^^" at the CERN LHC. 



D. Lepton angular distribution in the laboratory frame 



In this section, we give numerical results for the angular distributions of (.~ with respect 
to the electron beam axis computed with complete spin correlations between production and 
decay. Unlike the invariant mass distribution, this lepton angular distribution is crucially 
dependent on the production-decay spin correlations. Moortgat-Pick and Fraas [|I^ have 
found in their detailed study that the effect of the spin correlations for the lepton angular 
distribution amounts up to 20% for lower energies and the shape of the lepton angular 
distribution is very sensitive to the mixing in the gaugino sector and to the value of the 
slepton mass. These points can be confirmed by comparing two results displayed in Figure 16. 
In the scenario SI (left figure) with large selectron masses, the lepton angular distribution 
is forward-backward symmetric and larger in the forward-backward direction. On the other 
hand, in the scenario S2 the angular distribution is forward-backward asymmetric, maximal 
near cos 61- = but suppressed in the forward-backward directions. Furthermore, the size 
of the lepton angular distribution and the forward-backward asymmetry depends rather 
strongly on the CP phases. 



E. Triple momentum product 

So far, we concentrate mainly on the CP-even production-decay correlated observables 
which depend on the CP phases only indirectly. However, the initial electron momentum 
and two easily-reconstructible final-state leptons allows us to construct a T-odd observable 

OT = Pe- {Pi' y<Pi+). (63) 

with I = e, fi. This T-observable Ot [0 can be finite if there exist non-vanishing CP- 
violating or CP-preserving complex phases in the amplitude for the correlated process. 
If we neglect the heavy particle widths, the T-odd observable can be utilized to directly 
measure the CP phases or to constrain them. 

Due to the general property of the T-odd observable, we know that it should be given 
a linear combination of the CP-odd quartic charges and D4. Note that in the scenario 
SI with large slepton masses, both of them are proportional to a small suppression factor 
(s^ — 1/4) as noticed in Eq. (42). Therefore, the T-odd observable as well as the electron 
EDM can not give any significant constraints on the CP phases in the scenario SI. On 
the other hand, since the slepton exchanges give a large contribution in the scenario S2, 
there might be a relatively large value of the T- (CP-)odd observable for no-trivial phases. 
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Furthermore, since the fundamental structure of the electron EDM determined by the CP 
phases will be very different from that of the T-odd observable, two CP-odd quantities 
will play a complementary role in constraining the CP phases. In order to make a concrete 
comparison of them, we explore the exclusion region by the T-odd observable on the CP 
phases S2 in the scenario S2. Since we cannot estimate the precise systematic uncertainties 
mainly related with the detection quality, we neglect the systematic uncertainties but take 
into account only the statistical errors. In this case, the boundaries of the excluded region 
of the CP phases $^ and $i at the N„-a level for a given integrated luminosity satisfy the 
relation 



/ 



= — — (64 

2 \{UT)\'^atot 



where (X) = / X ^^^d^ / atot over the total phase space volume $. In determining the exclu- 
sion area we take into account two possible combinations of two final-state leptons;(e~, e"*") 
and {iJ,~ , n'^) , which is responsible for the factor 2 in the denominator. 

Figure 17 exhibits the excluded area of the CP phases by the electron EDM measurements 
at 95% confidence level (shaded region) and by the T-odd observable Ot (hatched region) 
at a 2-(T level with an integrated luminosity of 200 fb~^. One can clearly notice their 
complementarity role played by two independent CP-odd quantities. The T-odd observable 
enables us to exclude all the range of $^ for most values of $i except for $i = 0, tt, 27r. In 
addition, they are complementary in the sense that the electron EDM is an indirect physical 
quantity determined at a very low energy, which does not need to observe SUSY particles but 
the T-odd observable Or is a direct observable to measure the SUSY CP phases exclusively, 
which however requires producing neutralinos directly. So, we conclude that through our 
investigations the T-odd observable can be a very efficient and complementary quantity in 
constraining or determining the CP phases if the lightest and second lightest neutrahnos are 
pair-produced and unless the sleptons are too heavy. 



VII. CONCLUSIONS 

In this paper, we have investigated the associated production of neutralinos e+e~ X1X2 
accompanied by the neutralino leptonic decay X2 ~^ X?^^^ i taking into account initial beam 
polarization and production-decay spin correlations in the minimal supersymmetric standard 
model with general CP phases but without generational mixing in the slepton sector. The 
stringent constraints by the electron EDM on the CP phases have been also included in the 
discussion of the effects of the CP phases. 

First of all, we have described possible flavor-preserving - selectron, chargino and neu- 
tralino - mixings in the MSSM with general CP phases without generational mixing and 
applied them to the evaluation of the electron EDM to investigate its dependence on the 
phases. As a result, we have identified two typical scenarios; one has large selectron masses 
of the order of 10 TeV and the other relatively light selectron masses of 200 GeV and a large 
higgisino mass parameter. The first scenario allows the full range for the CP phases <1>^ and 
$^ relevant to the neutralino process while the second scenario allows a finite space for the 
CP phases. Employing the allowed space as the platform for further investigations of the 
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neutralino processes, we have obtained several interesting results, which can be summarized 
as follows: 

• The production cross section and the branching fractions of the leptonic neutralino 
decays arc very sensitive to the CP phases. As a result, the total cross section is very 
sensitive to the CP phases. 

• If the electron/positron masses are neglected, the initial longitudinal and transverse 

polarizations of the initial electron and positron beams lead to three CP-even distribu- 
tions and one CP-odd distribution, which can be studied independently of the details 
of the neutrahno decays. While they are sensitive to the real SUSY parameters, those 
observables, especially, the T- (CP-)odd observable Vn is (almost) insensitive to the 
CP phases in both scenarios. 

• The production-decay spin correlations lead to several CP-even observables among 
them we have studied the two-lepton invariant mass distribution, the lepton angular 
distribution, and one interesting T-odd (CP-odd) triple product of the initial electron 
momentum and two final lepton momenta. On the whole, we have found that the 
distributions are sensitive to the CP phases in the scenario S2 with relatively light 
selectrons and large gaugino compositions of the neutralinos. 

• We have presented the exclusion region of the CP phases $^ and $i by the T-odd (CP- 
odd) observable with the assumed integrated luminosity of 200 fb~^ at 2-a level. In 
comparison with the constraints from the electron EDM measurements, the constraints 
from the T-odd observable is complementary in that it constrains very strongly the 
phase $1. 

To conclude, the associated production of neutralinos e^e~ xtXi followed by the 
leptonic decays X2 ~^ Xi-^"^-^" is expected to be one of the cleanest SUSY processes 
to allow for a detailed investigation of the physics due to the CP phases in the MSSM. 
Therefore, if the neutralinos are produced at future e+e" colliders, the colliders will make 
it possible to measure or constrain the SUSY parameters and CP phases and so provide 
a complementary check for the existence of CP violation in the MSSM in the neutralino 
sector. 
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FIG. 2. (a) the allowed range for the phase vresus the higgsino mass parameter and (b) 
the allowed region of the phases and against the electron EDM constraints. The trilinear 
parameter \Ae\ is taken to be 1 TeV and it phase is scanned over the full allowed range. 
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FIG. 4. Neutralino mass spectrum; (a) m^o and (b) ni^o on the {^^, ^i} plane in the scenarios 



51 (upper part) and 52 (lower part). 
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(a) Scenario 1 



(b) Scenario 2 




FIG. 5. Gaugino compositions of the two lightest neutrahno states Xi X2 with respect to 
the CP phase with the phase $1 scanned over its full range in (a) in 51 and (b) 52. The open 
circles are for Xi and the filled diamonds for X2- 




FIG. 6. Dependence of the differential production cross section of the process e~^e~ xixi on 
the scattering angle Q in the scenarios (a) 51 and (b) 52 for five combinations of the values of two 
CP phases and <I>i. 
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(a) Scenario 1 
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(b) Scenario 2 
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FIG. 7. Dependence of the production cross section a{e~^e~ — > X2X1) on the {$^,<I>i} plane 
in (a) the scenario SI and (b) the scenario S2 for five combinations of the values of two CP phases 
$u and $1. 
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FIG. 9. Dependence of the left-right asymmetry Alr and the transverse distribution Pt on 
the scattering angle 9 in (a) the scenario SI and (b) S2 for five combinations of the values of two 
CP phases and $1. The upper two figures are for Alr and the lower two figures for Pt- 
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Normal Polarization in Scenario 2 



0.3 



0.2 



0.1 



0.0 



-0.1 



- (0. 0) 
(0,11/2) 

- (11/2,0) 
■ (1/2,11) 

(1,1/2) 



-1.0 -0.5 



0.0 

COS0 



0.5 



1.0 



FIG. 10. Dependence of the X2 normal polarization VjJ- on the scattering angle @ in the scenario 
S2 with relatively Hght selectron masses and a large value of = 700 GeV for five combinations 
of the values of the CP phases and $i. 
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FIG. 12. Contours of the branching fraction B{x2 ~^ Xi^^ 
plane in (a) the scenarios (a) SI and (b) the scenario S2. 
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(a) Scenario 1 
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(b) Scenario 2 




FIG. 14. Contours of the total cross section of the associated production of neutralinos 
e+e~ — X2X1 followed by the sub-sequential decay X2 ~^ Xi^^^~ on the plane of two phases 
and $1 in the scenarios (a) 51 and (b) S2. Here, i is either e or /x. 




FIG. 15. Two-lcpton invariant mass distributions of the correlated production process 
e+e~ X1X2 ~^ Xi'^~^^~ ill the scenarios (a) 51 and (b) 52 for five combinations of the val- 
ues of two CP phases and $1. 
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FIG. 16. Lepton angular distributions of the correlated process e'^e —^ X1X2 ~^ Xi(Xi^^^ ) 
in the scenarios (a) 51 and (b) 52 for five combinations of the values of two CP phases $^ and <&i. 
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FIG. 17. Excluded region of the phases {^f^, $1} by the electron EDM constraints (shadowed 
region) and by the triple product measurements with the integrated luminosity of 200 fb~^ at the 
2-a level. 
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